Model of Inhomogeneous Impurity Distribution in Fermi Superfluids 
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The standard treatment of impurities in metals assumes a homogeneous distribution of impurities. 
In this paper we study distributions that are inhomogeneous. We discuss in detail the "isotropic 
inhomogeneous scattering model" which takes into account the spatially varying scattering on the 
scale of the superfluid coherence length. On a large scale the model reduces to a homogeneous 
medium with renormalized parameter values. We apply the model to superfluid ^He, where porous 
aerogel acts as the impurity. We calculate the transition temperature Tc, the order parameter, and 
the superfluid density. Both A- and B-like phases are considered. Two different types of behavior 
are identified for the temperature dependence of the order parameter. We compare the calculations 
with experiments on ^He in aerogel. We find that most of the differences between experiments and 
the homogeneous theory can be explained by the inhomogeneous model. All our calculations are 
based on the quasiclassical theory of Fermi liquids. The parameters of this theory for superfluid 
''He in aerogel are discussed. 

PACS numbers: 67.57.Pq, 74.20.Fg 
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^\ . The standard treatment of impurities in a metal is based on averaging over the locations of the impuritiesi This 
' means that the probabiUty of a quasiparticle being scattered is independent of its location. To be definite, we call 
this the homogeneous scattering model (HSM). The pm'pose of this paper is to study the case where the impurity 
' distribution varies in space. We define a model of inhomogeneous scattering, which is as simple as possible and 
represents a medium that is uniform and isotropic on a large scale. We call it the isotropic inhomogeneous scattering 
model (IISM).^ The model is motivated by experiments on superfluid '^He in a porous aerogel, and all our numerical 
computations concentrate on this case. However, the model is independent of the pairing symmetry and therefore 
can equally be applied, for example, to s- or d-wave superconductors. 
' Liquid '^He is an ideal example of unconventional superfluid because it is naturally pure, it has a spherical Fermi 
surface, and its p-wave pairing state is well understood. Therefore it is of interest to study impurity effects in this 
superfluid. The addition of impurities to helium can be done by using porous aerogel so that typically 98% of the 
volume is occupied by '^He. It was found experimentally that the superfluid transition temperature in aerogel is 
reduced but remains sharp. Also other superfluid properties such as the superfluid density and the NMR shifts were 
M found to be modifled by aerogel. Many experiments studying this system have been made during the last seven yearsi^ 
The HSM is an attractive model for ^He in aerogel because of its simplicity. Essentially all theoretical calculations 
^ for this system are based on it,2iSiSiL2i24i2iii4iSiiiiiiiS4i& However, already the first comparisons to experiment showed 
O ■ that the HSM is insufficient quantitatively. In particular, it was found that the order parameter is more suppressed 
^ ' than the superfluid transition temperature Tc. A natural explanation for this comes from the fact that the scattering 
^ , in aerogel is not homogeneous, as has already been demonstrated using the IISMi Similar results have been reached 
r ""j ■ in calculations on unconventional superconductorsiii 

' Preliminary results of the HSM have been reported in Refs. andlT^ In this paper we present the IISM in detail. 
Our studies are based on the quasiclassical theory of Fermi liquids. We discuss the assumptions of this theory (Sec. 
" " and how these are satisfied for ^He in aerogel (Sec. IIIl)l . The inhomogeneous scattering model is introduced in 

Sec. IIVI For the case of p-wave pairing we introduce the order parameters of A- and B-type phases (Sec. We 
calculate several quantities including the critical temperature, order parameter, and superfluid density (Sec. lVl|l . The 
results are compared with experiments on ^He in aerogel (Sec. IVlfl . The equations of quasiclassical theory and the 
details of calculations are discussed in the Appendix. 



II. QUASICLASSICAL THEORY 



On a microscopic scale a pure system is described in terms of particles (conduction electrons or ^He atoms) and 
their interaction. Because the interactions are strong, this leads to a complicated many-body problem. 



The characteristic length in superconductivity or superfluidity is the coherence length. We define this quantity as 



iTTkB-LcO 

where v-p is the Fermi velocity. To be precise, we have used the superfluid transition temperature Tco of a pure system. 
The coherence length is typically much larger than the Fermi wave length Ap = 27r//cF- In '^He depends on pressure 
and changes from 16 nm at the solidiflcation pressure to 77 nm at the vapor pressure whereas Ap ~ 0.7 nm. 

The theory that is designed to work on the scale ^ '^f is the quasiclassical theory.i^i It treats the system as a 
dilute gas of weakly interacting quasiparticles. In quasiclassical theory all the many-body physics that takes place on 
the microscopic scale Ap is eliminated. It only appears through phenomenological parameters like the Fermi-surface, 
Landau Fermi-liquid parameters and transition temperature Tcq. 

Let us consider any external objects in the system. These objects are characterized by a strong potential, on 
the order of the Fermi energy. As a consequence the state of the system is modified in the vicinity of the object. A 
theoretical analysis of these atomic scale changes is again difficult because of the strong interactions between particles. 

In the quasiclassical theory the effect of external objects is twofoldiSi First, the phenomenological parameters 
discussed above are changed. These parameters are determined by processes on the Fermi-energy scale and are 
therefore of short range. Assuming that the surface area of the external objects times the atomic length scale (Ap) is 
a small fraction of the total volume, this effect is small and is neglected in the following. Second, the objects affect the 
low energy processes (energy ^ fcsTco) directly via a scattering of quasiparticles. The range of this effect is long, on 
the order of the coherence length [Eq. ^]. Therefore it leads to substantial modification of the superfluid properties. 

An important length characterizing the scattering is the mean free path £. The simplest case is to consider the limit 
Xp/£ 0. This is the quasiclassical limit where the Fermi wave length effectively disappears from the theory. First 
order corrections in Ap /£ lead to effects like weak localization, which are neglected here. In quasiclassical theory, the 
scattering is represented by a collection of "scattering centers." The main assumption is that the quantum interference 
between two scattering centers is neglected. Technically this can be achieved considering an ensemble average where 
the locations of the scattering centers are uncertain by a distance on the order of Ap or more. The size of one 
scattering center is limited by the condition that it has to be localized on the scale of ^q. The scattering properties of 
a center can be parametrized by scattering phase shifts S^^\ which are taken at the Fermi energy in the normal state. 
(For simplicity we label the different partial waves by a single index Z, but there is no need to restrict to spherically 
symmetric scattering centers. The scattering could also be spin-dependent, but it is also neglected here for simplicity.) 
Thus a complete description of the scattering needed in the quasiclassical theory consists of distribution functions 
ni{r) and scattering phase shifts sf^ of the scattering centers i — 1,2, . . .. 

Within each scattering center an exact quantum treatment is allowed in principle. However, because of uncertainty 
about the microscopic processes at surfaces, the phase shifts Sf^ cannot be calculated from flrst principles. Instead, 
one has to use some models, for example hard spheres]^ In some quantities the phase shifts only appear in certain 
combinations, for example, the transport cross section 

A ^ 

a = ^Y.(l + l)sin\6(^^'^~S^'^ (2) 
''f 1=0 

Let us try to clarify some consequences of the assumptions made above. The energy one usually is trying to calculate 
is on the order of /condCo- Here /cond ~ fceTc/Ap^o is the superfluid condensation energy per volume, and a typical 
volume ~ ^Q. Above the individual scattering centers were required to be small, cr ^ ^q. This implies that the energy 
/condfCo associated with a single impurity^S, is small in comparison. The typical number of impurities in volume is 
N ^q/(7, which for random impurities implies a fluctuation SN ^ \/N . The corresponding fluctuation in energy is 
by factor \/a/(,o smaller than /cond^o- This has to be neglected since there are other neglected contributions that are 
on the same order of magnitude. Thus impurity averaging in the quasiclassical approximation implies a scattering 
medium where the fluctuations in the impurity density are neglected. 



III. AEROGEL 

The structure of aerogel, as relevant for '^He experiments, is discussed in Ref. l2,Sl Here we repeat some main points. 
Aerogel consists of small Si02 particles of diameter ^ 3 nm, which are coalesced together to form a self-supporting 
structure. Experiments with '^He typically use aerogels with open volume fraction 98% or more. According to small- 
angle x-ray scattering measurements, there is a "fractal" range in the particle cluster up to a "correlation length" 
fa ^ 100 nm. Above this scale the structure looks homogeneous. Computer simulations give a picture of widely 
spaced aerogel strands implying a long mean free path > 100 nm. 



It seems reasonable that the quasiclassical description above can be apphed to hquid '^He in aerogeL The atomic 
layer on the Si02 surfaces occupies only one per cent of the volume and can be neglected for many purposes. (Magnetic 
properties make an exception because the susceptibility in this layer is much larger than in pure liquidi^) Since a major 
part of the liquid is within a coherence length from Si02r^ the scattering effect leads to a substantial modification of 
the superfluid properties. It also seems that the scattering from aerogel can be represented by incoherent scattering 
centers whose size is small compared to because random variations on the order of Ap are likely to develop already 
at much smaller distances. 

The smallest reasonable choice for a scattering center is a single Si02 particle (diameter 2Ra ^ 3 nm) . This is large 
in comparison to Ap: kpRa ^ 10. According to hard sphere phase shifts this means that only 1% of the scattering 
takes place in the s-wave channel, and 99% is left to higher partial waves in Eq. (jSJ. This dominant contribution of 
higher partial waves has several important consequences. 

Firstly, the phase shifts (S|''' are random numbers. This is because even the phase shifts (modulo tt) of a hard sphere 
with a fixed Ra are pseudorandom numbers for k^Ra S> 1. Adding to this the varying particle size, the surface 
roughness and varying orientations of touching neighbor particles, it is simply impossible that the result would be 
anything else but random. (The randomness, of course, is valid only for partial waves I that contribute essentially to 
scattering, i.e., for I < kpRa-) 

A consequence of the random phase shifts is that only the number of scattering phase shifts, or equivalently, the 
cross section a [Eq. |(5J)] is important in describing a scattering center. Thus a sufficient description of the scattering 
is obtained by specifying only cr^ and ni(r). 

A second important consequence of the large particle size is that both ai and ni(r) are independent of pressure. 
This is crucial for comparison with experiments, because when fitting is needed, it can be done at one pressure only, 
and the predictions of the model get fixed at all pressures. We argue as follows. It is reasonable to assume that the 
aerogel is independent of the hydrostatic pressure, implying that ni{r) is also. The Fermi wave vector changes by 
10% over the pressure range from the vapor pressure to the solidification pressure. This could induce some pressure 
dependence in the cross section ai. For example, hard-sphere a depends essentially on fcp in the region k^Ra '--^ li^ 
However, in the physical region k-pRa ^ 10 the dependence of a on fcp is very weak. Thus we conclude that both ai 
and ni{v) are independent of pressure. 

The large k-pRa is potentially bad news for theory because calculations that take higher partial waves into account 
are very complicated, see Ref. 0. The promising conclusion of these calculations is that, at least in some cases, 
the results including higher partial waves are almost identical to those including only s-wave under the following 
conditions: (i) one uses the same transport mean free path € and (ii) one uses either random phase shift J^^-* or fixed 
sm (5(0) w ^ in the s-wave calculation. Here we assume that this correspondence holds more generally. Thus we 
calculate only s-wave and present results for sin^ (5(°) — ^ . Finally, the randomness of the phase shifts also simplifies 
the numerical calculations since it implies that some components of the propagator vanish (see the Appendix). 

The use of several impurity species is important in studying anisotropic scattering where the preferred direction 
varies in space. In this paper we limit to isotropic scattering. In order to simplify the notation we select all scattering 
centers to have equal cross section, so that the only scattering parameters are a and the total impurity density n(r). 

IV. SETTING THE MODEL 

The simplest possible impurity profile n(r) is a constant. It implies a location independent (transport) mean free 
path I ~ [na)^^ . This homogenous scattering model (HSM) has successfully been used to model the impurities in 
superconductors^*^ For impure p-wave superfiuid ('^He) it has been used to calculate the critical temperature?^ 
the order parameter and superfiuid densitji^'^'^^, properties near the superfluid transitions^ properties in magnetic 
field fS'ifi'iSiiii& density of stateSfiSiii thermal conductivity-'^-^ and estim ate s tr ong-coupling eff ectsi^ 

The predictions of the HSM are compared to experiments in Refs. I37I9I I11II12|IiMTrI|23II28I andl29l see also below. 
Compared to pure '^He, the HSM gives the right tendency and can work even quantitatively in some cases, but there 
can be differences up to a factor of 5 in the suppression factors (see below). In some papers the experiments are 
compared with data that is calculated in the unitary limit of s-wave scattering, which seems to give better agreement 
than other phase shifts. As explained above, we believe this is misguided, and random or intermediate phase shifts 
should be used instead. The observed differences can more naturally be explained by inhomogeneous scattering as 
will be shown below, at least for some quantities. 

Real aerogel has voids where scattering is negligible. This can be modelled by an impurity density n(r) which 
depends on the location r. It is in principle possible to use a realistic n(r) for aerogel. This has the drawback that 
the computational volume should be large in order to get a representative sample, and this implies heavy numerical 
effort. Here we prefer the opposite limit of a simple model ?i(r). The simplest possibility would be a plane-wave 
variation n(r) = no + ni cos(q • r). A stronger version of this would be equally spaced scattering planes. In the limit 



FIG. 1: The sphere used in the IISM to approximate a unit cell of a periodic lattice. The drawing illustrates the boundary 
condition that a quasiparticle exiting the sphere is returned back at the diametrically opposite point. The shading depicts the 
spherically symmetric impurity density n{r). 



of very strong scattering in the planes this leads to isolated slabs The problem with all these models is that they 
are anisotropic. For example, the superfluid density would depend on the direction of the superfluid velocity Vs. In 
comparison to experiments one should use some average over the directions, but this neglects the process how the 
averaging really takes place by a nonuniform current distribution. 

The purpose of the isotropic inhomogeneous scattering model (IISM) (Ref. 0) is to incorporate a non-constant n(r) 
with spherical symmetry. We take a spherical volume of radius R and use an impurity density n{r) that depends only 
on the radial coordinate r. We assume that these spheres fill all the space. This last point is not strictly possible, 
but represents an approximation that is similar to using spherical approximation for a Wigner-Seitz unit cell. — The 
calculation in a single unit cell can represent states where the superfluid order parameter A{r) has the Bloch form 



where A^^^ (r) is a strictly periodic order parameter. In the present case the wave vector q is the imposed phase 
gradient that is related to the superfluid velocity Vg = (/i/2m)q deflned on a scale larger than R. 

At the surface of the IISM sphere we impose the boundary condition that an exiting quasiparticle is effectively 
returned to the sphere at the diametrically opposite point, see Fig. ^ In the case of current-carrying states, there 
has to be a phase shift which corresponds to the phase factor in Eq. Q. Otherwise the state of the quasiparticle 
(momentum, spin) is unchanged. The boundary condition can be expressed mathematically for the Green's function 
as 



see the Appendix for notation. The model reduces to the HSM in the limit that n{r) is independent of r. 

Because of the spherical approximation, it is worth while to properly define all large-scale (^ R) quantities. We 
deflne macroscopic quantities as averages of the corresponding local quantities over the sphere. 



In particular, corresponding to the local mass current density of the superfluid, j^, we have the macroscopic (js). The 
superfluid density )os is then defined by (js) = 'ps^s +0{v'^) with Vg = {h/2m)q as defined above. The average mean 
free path £ave is defined by = a{n{r)). 

The parameters specifying the scattering are the radius of the sphere R, the average mean free path ^avc, and the 
shape of the impurity density n(r). For n{r) we use two different analytic forms 



2(r) =2(°)(r)exp(iq-r), 



(3) 



g{k, Rr, e„) = exp (iq • Rtts) g{k, -Rr, e„) exp (-iq • Rtts) , 



(4) 




(5) 




(6) 



(7) 



with parameters j and b. Here the prefactors c and c' are determined by the parameters ^avo and cr (c, c' > 0). The 
functions in Eqs. lO and iQ are shown in the inset of Fig. |21 We call them void and cluster profiles, respectively, 
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FIG. 2: The B-phase order parameter [Eq. 0] in the IISM for Vs = at the temperature O.STco. The sohd hnes are for steep 
void impurity profile [Eq. Q with j — S], R = 2£ave, and ^o/^avo ~ 0.2. The dashed lines are for gentle void profile [Eq. @ 
with j = 2], R = £a.ve, and ^o/^ave ~ 0.1. The dash-dotted lines are for cluster profile [Eq. Q with j — 3, b — 0.1], R = ^ave, 
and ^o/^avc = 0.1. The inset shows the corresponding impurity profiles. 



because the former has strongest scattering at the boundary and the latter at the center of the sphere. Both profiles 
have zero derivative at r = i? in order to have a smooth impurity density everywhere. 

An attractive feature of the IISM is that the symmetries of a homogeneous system are preserved on a large scale 
3> R. Thus one can apply phcnomenological large-scale theories such as Ginzburg-Landau and hydrodynamic theories 
as for a homogeneous mediumiiSi*^ The only change is that the parameters of these theories are modified by the 
inhomogeneity. Some of these parameters are calculated below. 



V. ORDER PARAMETER 



We apply the IISM for superfluid ^He. The main assumption in addition to those already mentioned is the use of 
the weak-coupling approximation. The dipole-dipole interaction is neglected because it is unimportant on the scale 
of a few which we study here. The equations and some details of the numerical implementation are discussed in 
the Appendix. 

Even with all simplifications, the computational effort in the IISM is quite substantial. For example, the computer 
code has five nested loops in addition to the one needed for the iteration of the order parameter, and there are several 
stages of initialization, interpolation, and data collection. 

Based on calculations with the HSM, no new superfluid phases of '^He are expected in the presence of scatteringiS 
However, the order parameters of the A and B phases are modified by inhomogeneous scattering. The general forms 
can be deduced using symmetry arguments. The B phase order parameter for v^. = has the form 

A{r, 0, 9) = e'^R [a, (r)rf + A„ (r) {99 + 0<^)] , (8) 

where spherical coordinates (r, 0, 9) are used. The phase x ^^id the rotation matrix R are arbitrary constants. The 
calculation determines the real- valued radial Ar{r) and angular Aa{r) functions. These functions are shown in Fig. [3 
for three different scattering profiles n(r). We see that the order parameter is inhomogeneous, and is mostly suppressed 
in regions where the scattering is strong. For the void profile [Eq. ©] the order parameter components have maxima 
Aniax = Ar(0) — Aa(0) at the center and minima at the surface. The absolute minimum Amin is given by Aa{R). 
We define two different averages of the spatially varying order parameter. The usual average is given by 



AL = ^(A^ + 2A^). (9) 



An "NMR average" is defined by 



A^MR = ^(2A^-f6A,A, + 7A2). (10) 



This form can be justified as foUows. The frequency shifts of nuclear magnetic resonance in superfluid '^He are 
determined by the dipole-dipole interaction energy^i 



/d = <?d[|TrA|2 + Tr(AM)]. (11) 

The coefficient ^d is a phenomenological parameter that depends on a cut-off energy, which is on the order of the 
Fermi energy.'^'* According to the principles of Sec. ^ we assume it is not changed by the impurity. We calculate the 
dipole-dipole energy [Eq. (|11|) ] for the B-phase order parameter © and find 

/d = 4gdA^MRCOSi9(l + 2cosi?), (12) 

where ■& is the rotation angle of R. This has exactly the same form as in pure ^He-B, the only change being that the 
energy gap A is replace by Anmr- Thus the NMR properties within IISM are the same as in bulk liquid except for 
this renormalization. 

We point out that the order parameter in Eqs. (P^-- (|17|l is defined using the off-diagonal part of the mean-field self 
energy (jA4|l . It should be noted that the order parameter is not simply related to the energy gap in the excitation 
spectrum. Only in pure, homogeneous superfluid the energy gap equals Aavo = Anmr- The excitation spectrum in 
the HSM has been studied in Refs.|2 andlTM. but there are no studies yet for the IISM. 

The A-phase order parameter is more complicated since it is anisotropic. We select the cylindrical coordinates 
(p, 4>, z) so that z is along the anisotropy axis 1. The order parameter for = can be written as 

0, z) = e'-^d [Ap(p, z)p + iA^(p, z)0 A,(p, z)z] , (13) 

where d is an arbitrary constant unit vector. The functions A^ are real but now they depend on two coordinates p 
and z. They satisfy symmetry relations Ap(p, —z) = ^p{p, z), Acf,{p, —z) = Aij,{p, z) and Az{p, —z) = — A^(p, z). For 
the A phase we define 



AL = -(A^ + A| + A^) (14) 

A^MR = \{Al + Al-2Al). (15) 

In pure homogeneous superfluid Aavc = Anmr equals the maximum energy gap in the A phase. The dipole-dipole 
energy [Eq. I|ll(l ] in the A phase (|13|l is given by 

/d = -2.9dA^MR(d-i)', (16) 

and the justiflcation for Anmr 1115(1 is completely analogous to the case of the B phase. 

In the case of a flnite superfluid velocity (taken to be in the z direction) we limit our calculations to the B phase 
where the order parameter takes the form 

A(p, 0, z) = e'^S Apppp + A^^^4) + A,^zz + Ap^pz + A, pzp . (17) 

Now Api are complex and depend on p and z. They satisfy symmetry relations Apz{p, —z) = —A*p^{p, z), Azp{p, —z) — 
— A*p(p, z) and Aii{p, —z) = A*j(p, z), where i = p,(j), z. The mass supercurrent jg (r) has non-zero p and z-components 
inside the sphere. One case is illustrated in Fig. O It can be seen that the current density is smallest in regions of 
strong scattering. Naturally, the current is conserved, V • js(r) = 0. In the void profile, the current has to go through 
the scattering regions at the cell boundary, whereas in the cluster profile (not shown) the transport current can flow 
past the scattering region situated at the center. Regardless of the proflle, the averaged current is parallel to Vg and 
independent of the direction. Thus the superfluid density tensor /5s reduces to a scalar p^. 



VI. RESULTS 



We start by studying the critical temperature T^. We calculate the quantity Tc/Tco, the critical temperature 
relative to the critical temperature in the absence of scattering. In the HSM this depends only on the parameter 
^o/^- The dependence turns out to be the same as calculated by Abrikosov and Gorkov for s-wave superconductors 
in the presence of magnetic impurities^ This result was generalized to the nonmagnetic p-wave case in Ref. In 
the IISM Tc/TcQ depends only on the ratio ^o/^avc and the impurity profile n{r). In particular, it is independent of 
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FIG. 3: The local superflow js(/5, 2) in the IISM unit cell. The parameters are the steep void impurity profile [Eq. with 
j = 8], R = £ave = 5^0, T = 0.6Tco, Fl = 10 (pressure 15.6 bar), and q^o = 0.001. 




FIG. 4: The transition temperature as a function of average scattering rate ^o/^ave. The dotted line is from the HSM. Solid 
lines are for the steep void impurity profile [Eq. ® with j = 8], and correspond to R/£s,ve = 5, 1, and 2, from left to right. The 
dashed lines are for the gentle void profile [Eq. @ with j — 2] and from left to right R/ia.ve = 1 and 2. The dash-dotted lines 
are for cluster impurity profile Eq. Q with j — 3 and R = i^vc- The lower curve is with b — 0.1 and the upper one with fo = 0. 



the phase (A or B) and the phase shifts S^^\ The relative Tc is plotted as a function of ^o/^avc in Fig. 0] for different 
scattering profiles. One sees that in aU cases the HSM gives the lowest Tc- This is natural since the inhomogeneity 
imphes existence of regions where the scattering is less than the average, and in these regions the order parameter 
nucleates at a higher temperature. For example, for the steep void profile with R/iavc = 2 we find Tc — O.STco at the 
average impurity where superfluidity in the HSM is completely suppressed (^o/^avc = 0.28). Generally one concludes 
that the larger the inhomogeneity (in amplitude and in length scale), the larger Tc is obtained at a given ^o/^ave- One 
notices that the cluster profile with no background scattering (6 = 0) differs qualitatively from the other curves in 
Fig. ^ This profile is exceptional because it has long quasiparticle trajectories where the scattering is negligible. 

Let us now turn to the amplitude A of the order parameter. In the HSM A/fceTc depends mostly on the relative 
temperature T/Tc but also on the scattering rate S,o/i, the phase shifts S^^\ and the phase (A or B)^ The dotted lines 
in Fig. represent four different values of £,o/£ for B-phase with sin'^ 

In the IISM the parameter ^0/^ is replaced by ^o/^avc and n(r). The order parameter depends on the location r ||SJ|. 
In order to describe its temperature dependence we use parameters A^ax and Amin, and Aavo ®- The temperature 
dependence of the three characteristic numbers is shown by solid lines in Fig. We see that there can be strong 
variation since A^j^^ and A^^^^^ are quite different. The average A^^^ is considerably below the HSM curve that gives 
the same Tc- Therefore A^^^, is more suppressed than the critical temperature Tc- 




FIG. 5: Temperature dependence of the squared order parameter in the B phase. The sohd Unes give the minimum (A^i,j), 
the maximum (A^ax)i ^nd the average value [A^^j,, Eq. Q]. They are calculated in IISM for void impurity profile [Eq. 
with R = 2^avc. The left panel uses the steep profile {j = 8) and ^o/^avc = 0.275 and the right the gentle profile [j = 2) and 
Co/^avc = 0.16. Both cases have the same transition temperature (Tc = O.SlTco). The three dashed lines in both panels are 
calculated with the HSM. From larger to smaller they correspond to the pure limit {£,o/t = 0), the HSM with the same Tc 
i^o/f- = 0.076) and HSM with the same average scattering rate as in the IISM ((,o/£ — 0.275 in the left panel and 0.16 in the 
right panel). 



The temperature dependence of A^^^, varies. Let us first consider the case in the right-hand panel of Fig. where 
the shape of A'^^^{T) is concave near Tc. This is typical for a small inhomogeneity ampUtude. In our case it also means 
that the average scattering is small since the minimum scattering vanishes. We see that at low temperatures A^^^, 
agrees very well with the HSM curve that is calculated with the same average scattering rate. At higher temperatures 
A^.^,g deviates from this since the true Tc is higher than the one given by the average scattering rate. 

In the case of the left hand panel of Fig. the temperature dependence of A^^^ is quite linear near Tc- In fact, 
the linear range is wider (relative to Tc) than in the pure case or HSM. This happens when the amplitude of the 
inhomogeneity is large. Large inhomogeneity means large average scattering, thus the HSM results based on the 
average scattering rate are more suppressed than in the case of concave Al^^iT)- In the particular case of Fig.^Ia), 
the HSM result is quite small since £^o/£ — 0.275 is near the critical value of complete suppression. In this case A^^c 
is nowhere nearly approximated by the HSM result. However, the linear A^^g(T) can also appear in cases where 
the HSM with average scattering rate is not that suppressed, and provides a good approximation for A^.^,^ at low 
temperatures. This takes place at small R/^Oi where the proximity coupling between different regions tends to average 
out the inhomogeneity. In this limit the range of scattering rates /^avc where A^^^ C^) is concave near Tc seems to 
vanish. 

The concave shape of A^^^ (T) in Fig. IS^b) can be understood so that the different regions of an inhomogeneous 
sample have transitions more or less independently of each other: At Tc only A(r = 0) = Amax starts to grow but 
Aa(r = R) = Aniin and A,.(i?) (not shown) remain negligible until they start to grow at a lower temperature. In 
spite of the inhomogeneity, the onset of superfluidity indicated by Amax is very sharp giving a well defined Tc. 

The superfluid density ps is conveniently expressed relative to the density of the liquid, Ps/ptot- In addition to 
parameters discussed above, this depends on the Fermi-liquid parameter Tf in both the HSM and IISM. (In the IISM 
also other parameters could contribute, but they are neglected here.) Since ^ 10 is quite big, it has a strong effect 
on the results. The HSM results for "p^ in both the A and B phases have been calculated in Refs. andflU 

In the IISM we calculate ps for the B phase. Some curves are plotted in Fig. El As above, we compare the IISM 
with the HSM calculated with the same average scattering rate /^ave • The results of both models are close to each 
other at low temperatures when the scattering rate is small. At higher temperatures ps in IISM is larger because of its 
higher Tc. Thus Ps is relatively more suppressed than Tc. At larger scattering rates, ps in the HSM becomes smaller 
and is completely suppressed at £,o/£ = 0.28, while ps in the IISM stays finite much beyond that. All these results are 
qualitatively similar as discussed above for A^vo- The discussion of concave or linear temperature dependence in the 
case of A^^g cannot directly be applied to ps since this behavior is largely masked by the effect of Tf . The superfiuid 
density ps is also different from A^vo because of current conservation, which forces quite different current pattern in 
the void and cluster impurity profiles. 




FIG. 6: Superfluid density in the B phase at Fi = 10 (corresponding to ^He at 1.56 MPa). The uppermost dashed line is for 
pure superfluid. The dash-dotted lines are for the HSM with ^o/£ = 0.05, 0.10, 0.15, 0.20, and 0.25 in order of decreasing Tc. 
The solid lines are the IISM results corresponding to ^o/^avc = 0.05, 0.10, 0.40 and R = £a,ve- On the left we have used the 
steep void scattering profile [Eq. © with j = 8] and on the right the cluster profile with j = 3 and b = 0.1. 



In order to reduce the effect of Ff, we define a bare superfluid density by 

Ps ^ (1 + ^P?) (Ps/Ptot) 
Ptot 1 + (ps/Ptot) 



(18) 



In the case of the HSM this is independent of and in the case of the IISM the dependence is small, only a few 
percent over the whole pressure range . . .3.4 MPa. We calculate at an intermediate pressure of 1.56 MPa where 
Ff = 10. 

The order parameter and superfluid density can be expressed compactly by defining suppression factors 



Here is the bare superfluid density and Aq the order parameter, both in pure superfluid. The parameter t is the 
temperature relative to the transition temperature. Additionally, we wish to eliminate the parameter £avo, which is 
not directly measurable. This is achieved in plotting the suppression factors as a function of (Tc/Tco)^. One such plot 
is shown in Fig. |3 

The important conclusion from Fig. [T] is that the suppression factors for the A and B phases are almost identical, 
although there is a clear difference in the values of Aavc- This generalizes the result found previously within the HSMn^ 
Also, the difference in S^^ between Aavo and Anmr is small. The differences increase with increasing inhomogeneity, 
so that the suppression factors for A^^g and A^^j^ differ by « 5% in our extreme case R — 2£avo- AH these differences 
are rather small and, in order to simplify the plots, we present below 5^2 only for the B phase Aavc- 

Suppression factors for both A^ and as functions of (Tc/Tco)^ are plotted in Fig. |S| The HSM results are shown 
for comparison. They are all above the diagonal. It is clearly visible that the effect of inhomogeneity is to bend these 
curves down. This means that p^ and A^ are more strongly suppressed than Tc. With a steep profile and a large 
R/£ave one obtains a strong suppression and a wide flat region in the suppression curve. Superfluid density is clearly 
more sensitive to the inhomogeneity than the order parameter. 



VII. COMPARISON WITH EXPERIMENTS 



A comparison of the calculated transition temperatures with experimental results is given in Fig. IHl The pressure- 
independent £avc is determined so that the theoretical curves and the experimental data intersect at the point where 




FIG. 7: Suppression factor for the squared order parameter in different cases: B pfiase A^vc (solid lines), B phase A^^r (dotted 
lines), A phase A^ve (dashed lines), and A phase A^mr (dash-dotted lines). The five sets of curves correspond to reduced 
temperatures t = 0.2, 0.5, 0.7, 0.8 and 0.9, from top to bottom. Other parameters are steep void profile [Eq. @ with j — 8] 
and R = £ave. The inset shows double magnification of the part indicated by the rectangle. 



Tc/Tco = 0.7. For the samples measured in Manchester [Fig. ^p] the HSM gives a fairly good fit but a better 
correspondence is obtained with IISM by using R = ^avc- For samples measured in Cornell and Northwestern [Fig.Eti] 
the closest curve is the IISM with R = 2£ave and j ^ 8 m Eq. This represents roughly the best fit obtained by 
the IISM. 

Particularly interesting are the samples A and C measured in Cornell [Fig.lHIa)], which also have been studied using 
small-angle x-ray scattering. The values measured for the aerogel correlation scale are 130 and 84 nmm^ These are 
slightly smaller but on the same order of magnitude as our closest curve values R = 210 and 140 nm, respectively. 
This supports the view that the better fit obtained by the IISM as compared to the HSM is not only due to more 
fitting parameters, but a more realistic modelling of the structure of the aerogel. 

Some experimental points for both the order parameter and superfluid density are shown in Fig.|Sl Each marker (V, 
A, * etc.) identifies a data set that corresponds to one aerogel sample at a fixed reduced temperature t. The different 
points in each data set are obtained from measurements at different pressures. The superfluid density is measured 
by torsional oscillator and the interpretation of the experiments is rather straightforward. The order parameter is 
measured by NMR. This is possible because the frequency shift depends on the dipole-dipole energy, Eqs. ((T^ or lfTC|l . 
For interpretation one has to know whether the phase is A or B. The Stanford point (A) is based on seeing similar 
suppression in both phasesJ^^ The Northwestern data is drawn here by assuming A phase and d and 1 perpendicular 
to magnetic field, as they are in bulk ^He-AiSii2& However, it is not clear if this original assumption is correcti^ In 
case of the B phase the analysis would be more complicated because it would require an analysis of the texturei2^'2L22 

One can notice that the experimental data for the order parameter and superfluid density are qualitatively similar. 
They all fall clearly below the diagonal, and appear to extrapolate to S* = already at a finite Tc/Tcq. This is in 
clear disagreement with the HSM results that are all above the diagonal. It can be seen that much better agreement 
is achieved with the IISM. Ideally each experimental data set should fall on some of the theoretical lines. As one can 
see, this is not quite the case in the IISM. A reasonable over-all fit to suppression factors and critical temperatures 
is achieved with the void impurity profile l|njl with j = 2 and R « l.S^avc, but the choice depends on the properties 
one wishes to emphasize. The agreement could be improved by allowing for a pressure-dependent mean free path £. 
However, the most obvious reason for the remaining differences is that the scattering profile in the IISM has only one 
length scale whereas real aerogel has a wide distribution of length scales. 

VIII. CONCLUSIONS 

We have presented the isotropic inhomogeneous scattering model. We claim it is the simplest model of inhomoge- 
neous scattering that is consistent with large-scale isotropy. Unfortunately, the computations needed are much more 
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FIG. 8: Suppression factors 1)19^ for the squared order parameter [Eq. 0] (left) and for the bare superfluid density l|18|l (right). 
The solid and dashed lines are for void impurity profile ^ with steep {j = 8) and gentle {j = 2) slopes, respectively. The 
dash-dotted lines are for the cluster profile Q with j = 3 and b = 0.1. The dotted lines are calculate using the HSM. For 
subplots a-d the temperatures from top to bottom are t = 0.2, 0.5, 0.7, 0.8, and 0.9 for the IISM. For HSM the upper line is 
for f — > and the lower one for t — > 1. The subplots e and / compare the different profiles at t — 0.8. The radius of the unit 
cell is mainly R = ^ave except for subplots c and d where the solid lines are for 7? = 2i'avc and for subplots e and / where the 
uppermost solid line is for R = 0.5^ave and the lowest solid and dashed lines are for R — 2^ave. Experimental points for the 
order parameter are from Refs. l2^ andl3^ with t = 0.9 (□) and t = 0.5 (+) and|^ (A, t — 0.65), see discussion in the main 
text. The measured superfluid densities are from Manchester (Ref. |2^ with t = 0.5 (x), t — 0.8 (o), and t = 0.9 (0) and from 
Cornell sample A (Ref. 153), with t = 0.5 (V), t = 0.8 (A), and t = 0.9 (*). 



demanding than in the homogeneous scattering model. The model itself is independent of the pairing symmetry and 
thus could be used to study the effect of inhomogeneous impurity distributions in superconductors. 

When applied to "^He in aerogel, the IISM gives better agreement with experiments than the HSM. We emphasize 
that this is not solely due to the IISM having more free parameters than the HSM. On the contrary, the inhomogeneity 
of aerogel is the most natural if not the only way to explain the differences between measurements and the HSM. 
The fitted parameters of IISM are in reasonable agreement with measurements and simulations on the structure of 
aerogel. The fit given by IISM is not perfect, though. The main problem with the IISM is that it contains only one 
length scale whereas real aerogel must have voids of various sizes. 
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FIG. 9; Comparison of the calculated transition temperatures with experiments. The solid and dashed lines are calculated 
using the void impurity profile Q with steep {j — 8) and eentle {j — 2) slopes, respectively. The dotted lines are from the 
HSM. In panel (a) the experimental points are from Refs. fj (•, Cornell sample A), 2^ (A, Cornell sample C), andl24l 

andlii (□). Here £ = 213 (314) nm for the HSM, R = 24vo = 140 (206) nm for j = 8 and 7? = 2£ave = 234 (342) nm for j = 2. 
The values in the brackets denote the upper curves. In panel (b) the measured points are from Ref.|2§ (♦, ▼) and we have used 
£ = 251 (324) nm for the HSM, R = ^ave = 132 (170) nm for j = 8, and R = ^ave = 195 (252) nm for j = 2. All experimental 
points are for aerogel with nominally 98% open volume except ♦, which are for 99%. Small-angle x-ray scattering gives the 
aerogel correlation lengths = 84 and 130 nm for samples denoted by ▲ and •, respectively (Ref. l2l?) . 



APPENDIX A 



Here we present the relevant equations of the quasiclassical theory and discuss their numerical solution. The central 
quantity is the quasiclassical propagator 5.^^ It is a 4 x 4 matrix, whose components can be represented using the 
Pauli spin matrices as 



(Al) 



In equilibrium the propagator g(k, r, e^) depends on the direction of momentum k, on the location r, and on the 
Matsubara energy — ■nk^T{2m — 1), where m is an integer. The propagator is determined by the Eilenberger 
equations, 



[iem-Ta - CT, g] + ihy^'k ■ VrQ = 0, 
99 = -1- 



(A2) 
(A3) 



Here [A, B] = AB — BA denotes a commutator and t; denote Pauli matrices in the Nambu space; = a^i^l. The self 
energy matrix a = a^ii + o^imp consists of mean field and impurity contributions. For spin-triplet pairing the former 
part has the form 



CTmf(k, r) 



!^(k, r) + i/(k, r) ■ a A(k, r) • ma;2 

iCT2 A* (k, r) • a i/(-k, r) - a.^v'i-'k, r) ■ aa^ 



(A4) 



The real- valued functions v and u depend on the Fermi liquid parameters Fj^'"' . Omitting F^, F^ and higher coefficients 
we have 1/ = and 



(A5) 



where (• • ■ denotes angular average. The off-diagonal part of CTmf defines the order parameter, which for p-wave 
pairing is A = A • k. It is determined by the weak-coupling self-consistency equation 



A(k,r)ln— +^fcBTV 



^^-3((k.k')f(k',r,e,„)^ 



= 0. 



(A6) 



The impurity self energy (Timp is related to the forward-scattering t matrix via 



a-imp(k, r, £„) = n(r)i(k, k, r, e„). 



(A7) 



where n(r) is the impurity density. The usual derivation of this formula assumes a uniform impurity density n = consti 
However, it is valid for an arbitrary function n(r) as long as it allows the impurity positions to be uncertain on the 
scale of the Fermi wave length Ap- We assume that only s-wave scattering is important. This allows us to write the 
t-matrix equation directly for the impurity self-energy 

) = n{r)vl + nNpv (^g{k, r, em)y^ ^imp{r, Cm)- (A8) 

Here 2iVF = m* kp / n'^ fi^ is the total density of states at the Fermi surface. The effective mass m* is related to the 
atomic mass 7713 by to* = (1 + ■|Ff)TO3. The scattering potential v is related to the scattering phase shift J^^-' by 
V ^ -(l/7r7VF)tan5(o). 

In the calculation we use the following symmetry of the propagator and self energy 

[u{k, r, e™)]^ = -T2u{k, r, em)T2- (A9) 

This limits the nonzero components of g to 10, and allows one to solve the Eilenberger equation with the "multiplication 
trick" of Ref. 40 using only five components. In the current-free case all equations are immediately compatible with 
the symmetry in Eq. (jA9|) . In the case of nonzero current a possible problem could arise from the fact that in iterating 
the impurity equation (|A8|) the product (5)^(9)^ might not be proportional to a unit matrix 1. The inconvenient 
terms always appear with odd powers of v. Considering that the phase shifts in aerogel are random, such terms 
average out and we can use the symmetry of Eq. I)A9|I also for nonzero current. 

By taking into account all the symmetries the impurity self energy in the IISM has the form 

/ N f 16 + iS ■ a R • 010,2 ^ , 1 t 

CimplI", Cm) = • T,* -r , -s ] + wl, (AlO) 

where R, 6, and 5 can be expressed analytically as functions of the propagator {g)^^ The term w is of no interest 
since it drops out in the Eilenberger equation l)A2|l . In cylindrical coordinates {p, (j), z) the only dependence appears 
in the vector R in the A phase, where R = Rq{p, z)e"*z and Rq is real. The terms 6 and S are always real. In the A 
phase ^ = and in the B phase d = (J^c/) and R = Rpp + RzZ. 

From the propagator g I|A1|I one can calculate the supercurrent density js using 

j,(r) = 2^TO3i;FiVFfcBr^(k(?(k,r,e^)^^. (All) 

The superfluid density ps is obtained by calculating the averaged current at small values of Vg . 

The first step in the calculation is to give initial values for the fields A, (limp, and ly. Normally we used values obtained 
from the HSM. The iteration of these fields was started by transforming the data from the cylindrical/spherical 
coordinates to "trajectory coordinates." These are Cartesian coordinates with one coordinate, u, along the trajectory 
and two others (6 and t) specifying the position of the trajectory in the unit cell. Next the Eilenberger equation was 
solved along the trajectories. This was done by first calculating an unphysical exponentially growing propagator along 
each trajectory using the fourth order Runge-Kutta method. The end values at the cell boundary were saved to be 
used as starting values for next step in iteration. An exponentially decreasing solution was deduced using symmetries. 
The bounded physical propagator was obtained as a commutator of these two unphysical solutions as explained in 
Ref. I40L The solution of the Eilenberger equation was repeated for all Matsubara frequencies \em\ < ^n- This was 
repeated for all trajectories. In the general case it represents loops for coordinates b, t and the angle a between the 
trajectory direction and the z-axis. For the B phase without a current only the impact parameter b with respect to 
the center of the sphere was actually needed because of symmetry. 

In the next step better estimates for A, CTimp, and v were calculated. For that the propagator was converted from 
trajectory coordinates back to cylindrical/spherical coordinates and the required angular averages were calculated. 
The contribution from higher Matsubara frequencies \e,n\ > ejv was approximated by a Ginzburg-Landau form- with 
Matsubara sums evaluated using the Euler-MacLaurin formula. Using the updated values of the fields, the process was 
started from the beginning. The boundary condition was used to calculate the value of the exponential propagator 
at the initial point of the trajectory from the value at the final point stored on the previous round. This loop was 
repeated until the fields converged. 

In the numerical algorithm the number of discretized points in cylindrical coordinates pi and Zj was approximately 
40 for both in the range from to R. The discretization step in the Cartesian trajectory coordinates was approximately 
the same. Simple interpolation formulas were used for the transformations between coordinate systems. The angle 
a was typically discretized by eight Gaussian points in the range < a < n/2. The number of positive Matsubara 



frequencies em used was typically less than 20 for temperatures above O.lTcO- We made also more accurate test 
calculations. The qualitative behavior remains the same, but there are inaccuracies on the order of 2% in the results 
presented here. 
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